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This paper discusses the progress of work to model high-speed supersonic reacting flow. The purpose 
of the work is to improve the state of the art of CFD capabilities for predicting the flow in high-speed 
propulsion systems, particularly combustor flowpaths. The program has several components including 
the development of advanced algorithms and models for simulating engine flowpaths as well as a funda- 
mental experimental and diagnostic development effort to support the formulation and validation of the 
mathematical models. The paper will provide details of current work on experiments that will provide 
data for the modeling efforts along with the associated nonintrusive diagnostics used to collect the data 
from the experimental flowfield. Simulation of a recent experiment to partially validate the accuracy of 
a combustion code is also described. 


1 Introduction 

The design and development of a propulsion system is accomplished with several levels of ana- 

lytic tools, ground-based testing and finally flight. This paper will discuss the progress of work 

to improve the state of the art of CFD capabilities for predicting high-speed reacting flow in en- 

gines. The research has several components including the development of advanced algorithms 

and models for simulating engine flowpaths as well as a fundamental experimental and diagnostic 
development effort to support the formulation and validation of the mathematical models. The 
paper will provide details of current work on experiments and the associated non-intrusive diag- 
nostics that will provide data for the modeling efforts. This effort will lead to the development 
of phenomenological models for Reynolds averaged Navier-Stokes codes, subgrid scale models 
for large-eddy simulation techniques, and reduced-kinetics models that can then be applied to the 
design of high-speed propulsion systems. 



2 Axisymmetric Coaxial Free Jet Experiments and Related Design Computations 

Computational fluid dynamics (CFD) methods that employ the Reynolds-averaged Navier-Stokes 
(RANS) equations are widely used in the design and analysis of hypersonic airbreathing engine 
flow paths. These methods require models of statistical quantities of the turbulence in their de- 
velopment which have to be empirically calibrated and validated. In particular, new models for 
turbulent Schmidt and Prandtl number, as well as for turbulence-chemistry interactions, are re- 
quired [1], While suitable data is available for low-speed flows with combustion, sufficient data 
is still lacking in the area of supersonic combustion. Goyne et al. report measurements using 
particle-imaging velocimetry of mean streamwise velocity in a dual-mode scramjet [2]. At the 
NASA Langley Research Center, several data sets have been acquired in a H 2 fueled supersonic 
combustor using the coherent anti-Stokes Raman spectroscopy (CARS) technique [3] and the dual- 
pump CARS technique [4,5]. The former technique gave temperature only while the latter gave 
both temperature and composition. Data included both mean flow and turbulent statistics, although 
the uncertainty in the latter was limited by both instrument precision and number of measurements 
from which the statistics were formed. International work in this area includes measurements in 
scramjet combustors conducted at ONERA (France) and DLR (Germany) using CARS [6], and 
other non-intrusive techniques. 

Available data are limited to only a subset of the important variables (which are temperature, com- 
position, density and velocity) in a limited number of geometries, and turbulence data are even 
rarer. To meet the need for more data, an Interferometric Rayleigh Scattering technique for mea- 
suring instantaneous velocity is being developed at the NASA Langley Research Center to com- 
plement the Dual-Pump CARS technique [7, 8]. This technique uses Rayleigh scattering from one 
of the CARS laser beams to measure velocity in the same instant that CARS measures temperature 
and composition. Details on the techniques are given in a later section. Analysis of streams of such 
instantaneous measurements allows formation of the statistical quantities (means, variances, and 
covariances) required by modelers. Experimental configurations to provide suitable flow fields are 
being developed. An axisymetric coaxial free jet was selected since it provides the good optical 
access required for the Rayleigh technique, and symmetry can be taken advantage of to minimize 
the number of required spatial points. In order to form accurate statistics, large numbers of mea- 
surements are needed at each given point. The use of symmetry allows the total number of mea- 
surements to be made manageable. Another requirement of the flow field is that it should provide 
data that is relevant to both H 2 - and hydrocarbon-fueled hypersonic air-breathing engines. The 
experiment should be capable of providing flows of various Mach numbers to establish the effects 
of compressibility. Finally, it is desirable to be able to create both supersonic combusting flows in 
which the flame is attached to the burner (flame held) and flows in which the flame is detached, 
since both types of flow can occur in hypersonic engine combustors. Figure 1 shows a schematic 
and photo of the coaxial free jet experiment. Figure 1(a) shows the burner and nozzle, sectioned 
along the axis, without bolts, tubes, spark plug and other fittings shown; Figure 1(b) is an image 
of the facility near the nozzle exit during operation, showing the supersonic jet of vitiated air and 
the laser beams of the CARS-Rayleigh system. The facility consists of a water-cooled combustion 
chamber (burner), a silicon carbide nozzle (sonic convergent, or supersonic convergent-divergent 
for M=1.6 and M=2), with an exit diameter of 10 mm and a coflow nozzle. Various combinations 
of gas flows to the burner are possible. In one set of experiments, H 2 or CH 4 fuel, air and 0 2 
are reacted to provide vitiated air at various temperatures (dependent on flow rates). In these ex- 
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Figure 1: The supersonic jet combustor and nozzle, (a) Sectional view (without bolts, tubes, etc.), 
(b) Image during CARS-Rayleigh optical system data acquisition. 
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Figure 2: OH contours showing the sensitivity of the flame location to the presence of angled co- 
flow injection. 
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Table 1 : The experimental test matrix 


periments, the coflow is of unheated H 2 or CH 4 , and the result (depending on temperature) is a 
mixing and reacting coaxial jet flow. In another set of experiments, the gas flows to the burner are 
II), O') or N- 2 , and air in such a ratio that the products contain excess II). The coflow is of air 
and the result again is a mixing and reacting coaxial jet flow. Depending on temperature and on 
the co-flow rate, the flame may be held at the annular base region formed between the central and 
coflow nozzle exits where the flow recirculates. Non-reacting cases can also be considered in the 
case where there is no excess H 2 in the center jet and no fuel in the co-flow. 

As part of the design effort, a parametric study of the proposed geometry and run conditions was 
conducted using CFD to determine the sensitivity of important experimental quantities, e.g. tem- 
perature, species, density and velocity, to a number of variables including: inner jet temperature, 
flame holding base size, co-axial jet injection angle, inner jet Mach number and sensitivity to com- 
putational values of the turbulent Prandtl and Schmidt numbers. An example calculation using the 
Vulcan CFD code [10] is shown in Figure 2, which shows the sensitivity of the flame location to 
the presence of a co-flow. This flow is with a Mach 2, hot hydrogen-rich inner jet without and with 
cold air co-flow. 

Once the experimental geometry was finalized, CFD was used to help define the experimental test 
matrix. This matrix, given in Table 1, includes variations in the inner jet Mach number, the tem- 
perature of the inner jet, the type of fuel (hydrogen or methane), the equivalence ratio of the fuel 
and air and the location of the fuel and air streams (inner jet or co-flow jet). The CFD calculations 
completed in this step provide a baseline to compare with improved models planned for develop- 
ment using data from this and other experiments in the program. Figure 3 contains tables and flow 
images showing the various types of flames that are observed in the actual experiment. Two types 
of images are shown: visible light (true color images) with a long exposure time, and false color 
infrared (IR) images acquired in the 8-9 micron (long wave) region at an exposure time of 10 
ms. The tables contain information pertaining to the state of the flame: no flame, detached flame, 
flame held at the base or at the external coflow boundary, and if the flame holding is marginal, i.e., 
at the point of extinction. Figure 3(a) pertains to cases with AT 2 -vitiated air in the center jet and 
II) in the coflow. The center jet sensible enthalpy is varied by adjusting the flow rates from that of 
Mach 5.5 flight to that of Mach 7, and the exit Mach number of the center jet is either 1.0 or 2.0. 
The coflow is subsonic at the nozzle exit and the equivalence ratio (the ratio of H 2 coflow rate to 
that required to consume all the 0 2 in the center jet, 0) is either 0.5 or 1. (It is not implied that the 
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Figure 3: Flame state for a matrix of supersonic operating conditions at sonic and Mach 2 pressure 
matched exit conditions. Mixed infrared (false color) and visible light images (two right-most im- 
ages). (a) Vitiated air center jet with subsonic hydrogen air coflow at an overall equivalence ratio <$>. 
( b) Hydrogen rich center jet with sonic coflow of air. 


coflow H -2 reacts only with the 0 2 in the center jet.) For M exit = 1, the flame is held at the exterior 
boundary in all cases, whereas, for M exit = 2, in some cases there is no flame and in others the 
flame is detached (stands off from the nozzle exit), but there are no cases with flame holding. With 
a detached flame, the trend with increasing center jet enthalpy is for the flame to move towards the 
nozzle exit. Figure 3(b) pertains to cases with excess H 2 in the center jet, either 33% or 50% by 
volume of the jet flow being H 2 . The exit static temperature is increased in increments of 100 K, 
and the exit Mach number is either 1.0 or 2.0. The coflow of air is sonic and pressure matched 
at the exit. In some cases, no flame is observed; in others, detached flames or flames attached at 
the base. The trend with increasing T exit and decreasing M exit is from no flame to detached flame 
to base-attached flame. Thus, a full range of flame behavior is attained which is ideal for model 
development and code validation. 

Work has just begun to simulate the axisymmetric coaxial free jet experiment utilizing both RANS 
and LES codes. This work has not yet been completed. Work has been completed, however, 
on a coaxial jet mixing experiment that preceded the current work, and this experiment has been 
simulated using RANS codes. Details of the experiment are described in two earlier papers [11,12]. 
Low-density helium, which serves as a simulant of hydrogen fuel, was chosen to allow detailed 
studies of mixing without chemical reaction. Oxygen is added to the helium jet as a diagnostic 
aid for an oxygen flow-tagging technique (RELIEF). Several methods are utilized to characterize 
the flow field including Schlieren visualization, pitot pressure, total temperature, and gas sampling 
probe surveying, and RELIEF velocimetry. A schematic of the coaxial jet configuration is shown 
in Figure 4. The rig consists of a 10 mm inner nozzle from which helium, mixed with 5 percent 
oxygen by volume, is injected at Mach 1.8 and an outer nozzle 60 mm in diameter from which 
coflowing air is introduced also at Mach 1.8. The velocity ratio between the two jets is 2.25, the 
convective Mach number is 0.7, and the jet exit pressures are matched to one atmosphere. The 






Static pressure tap 


- Center jet 
pressure tap 

P ref, CJ 


Helium / 5% Oxyger 
or Air T t r . 


Plenum 
pressure tap 
P ref, coflow 


Air 

r t, coflow 


All dimensions in mm 


Figure 4: Coaxial jet assembly cross-section 


Figure 5: Schlieren image of coaxial jet mixing (conical extension cap removed) 
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Figure 6: Experimental setup of the combined CARS - IRS system 


resulting flow downstream of the nozzles can be seen in Figure 5, which shows a Schlieren image 
of the flowfield. The development of the mixing layer between the central helium jet and the 
air jet can be seen along with the shear layer development between the air jet and the surrounding 
quiescent laboratory air into which the air jet exhausts. Shock-expansion wave structure emanating 
outward from the centerbody nozzle lip can also be seen. Inward propagating waves from the inner 
lip, due to the finite thickness of the lip (0.25 mm), can be observed in the air jet once they pass 
through the helium jet. These waves are not visible in the helium jet due to the low refractive index 
within the center jet. A third wave can also be observed emanating inward from the outer nozzle 
lip and traversing both the air jet and the helium jet. The detailed data taken from this experiment 
along with the imaged shock-expansion structure provide an excellent case for assessing the ability 
of a reacting flow code to simulate complex mixing processes in a high-speed flow. 

3 Simultaneous CARS and Interferometric Rayleigh Scattering 

It was recently reported [13] and is summarized here, the combination of a dual-pump coherent 
anti-Stokes Raman scattering system with an interferometric Rayleigh scattering system (CARS- 
IRS) to provide time-resolved simultaneous measurement of multiple properties in combustion 
flows. Time-resolved simultaneous measurement of temperature, absolute mole fractions of N 2 , 
0 2 , and II 2 , and two components of velocity in a Hencken burner flame were performed to demon- 
strate the technique. The experimental arrangement of the combined system is shown in Fig. 6. For 
the measurement of temperature and the absolute mole fractions of N 2 , 0 2 , and II 2 a dual-pump 


CARS method [4] is used. The system uses spectrally-narrow green (injection seeded Nd:YAG at 
532 nm) and yellow (NB Dye laser at 552.9 nm) laser pump beams and one spectrally broad red 
laser (BB Dye laser at 607 nm) beam as the Stokes beam. 

The beams are combined at the focusing point of a spherical lens Lc (focal length of 410 mm) in 
a folded BOXCARS geometry to probe Raman transitions of N 2 , 0 2 and H 2 . The input beams 
plus the coherent blue signal beam at 491 nm are collected and collimated by another spherical 
lens with the same focal length as Lc. The three input beams are captured in a beam dump while 
the blue signal beam is passed to a spectrometer. The CARS signal, which is a spectrally broad 
blue beam that contains N 2 , ()■> and II> spectra, is analyzed by a spectrometer and recorded by 
the CCD1 camera. The shape of these spectra provides information on the temperature while the 
relative intensities of these spectra provide a measure of the relative mole fractions. The spectrum 
is fit with a theoretical model to determine the temperature and mole fractions. 

The velocity measurement is performed simultaneously using an interferometric Rayleigh scatter- 
ing measurement system [7]. The same pulsed, seeded green laser beam employed in the CARS 
technique is used as a narrow-band light source for the Rayleigh scattering system. The receiving 
optics for Rayleigh scattering are designed to capture the Rayleigh scattered light from the green 
pump beam in the measurement volume while preserving the scattering angle information, and to 
mix it together with the unshifted light of the laser before it is passed through the interferometer. 
In this way more than one component of velocity at multiple points can be measured in one in- 
terferogram [7, 14]. The velocity components measured are those that bisect the angles of laser’s 
incidence and the collection optics. Since two directions are collected, two velocity components 
are measured. The Rayleigh scattered light from the two measurement directions are combined 
with a polarizing beam splitter (PBS) and directed through a Fabry-Perot etalon and onto a CCD 
camera for analysis (CCD2). The particular setup used here gives a range of measurable velocities 
from 0 to ~3 km/sec up to temperatures of about 2500 K. The CARS spectra and the Rayleigh 
spectra are acquired simultaneously by synchronizing the cameras with the green laser Q-switch 
pulse at 20 Hz. The spectra are subsequently processed, as described in reference [13] to determine 
the temperature, composition and velocity. To demonstrate the method, experiments were carried 
out in a Hencken burner, which provides an adiabatic H 2 -a ir atmospheric pressure flame. The 
flame was stabilized with a co-flow of N 2 . This flame was used because it provided a challenging 
high-temperature measurement environment while producing a known (near-zero) velocity. 

Simultaneous CARS and Rayleigh scattering spectra up to 1610 K are shown in Fig. 7 (a-c) and in 
Fig. 7 (d-f), respectively: (a) and (d) are in the co-flow of N 2 , (b) and (e) are in a high temperature 
flow containing N 2 and 0 2 , and (c) and (f) are in a high temperature flow containing N 2 and a low 
proportion of 0- 2 . The CARS spectral plots show the experimental data, the theoretical fits, and 
also the residual between them. These spectra were used to calculate the rotational-vibrational gas 
temperature and the mole fraction of H 2 , 0 2 , and N 2 . 

For the Fabry-Perot spectra shown in Fig. 7, only one fringe order from each spectrum has been 
analyzed for this paper. The narrow shaped peak, more visible in the spectra of higher temperature, 
is the reference laser frequency (no Doppler shift). The broader spectrum is the Rayleigh scattered 
light. Each figure contains two Rayleigh and two reference spectra corresponding to the two col- 
lecting directions V 1 (black trace) and V2 (red trace of smaller amplitude), differing by 34 degrees. 
Two-component velocity measurements in the range of 300-1600 K showed near- zero velocities 
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Figure 7: Simultaneous spectra of CARS (left) and Rayleigh scattering from two viewing directions 
(right). The CARS spectra are both data and fits of theory to the data. The Rayleigh spectra are data 
points connected by lines 


(< 30 m/s) and standard deviations of 30-40 m/s. These errors are about 1% of the dynamic range 
of this measurement system (3,000 m/s). Measurement of one velocity component was shown to 
be possible at up to ~2,400 K with this system. These accuracies and precisions are within the 
desired range required for the planned supersonic combustion experiments where velocities will be 
up to 1,500 m/s. However, it is anticipated that future improvements in both hardware and software 
will allow these errors to be reduced by a factor of two or more. 


4 Simulation of Mixing and Reacting Flows 

Improved numerical simulation of high-speed propulsion systems and engine performance relies 
on the development of enhanced codes with an increased capability to model turbulence, turbulent 
mixing, and kinetics. A significant amount of the design effort utilizes Reynolds averaged Navier- 
Stokes (RANS) codes for flowpath prediction. RANS design codes utilize phenomenological mod- 
els to describe the turbulence field, fuel-air mixing, kinetics, and the interaction of turbulence and 
kinetics. There is an increased reliance on the use of large eddy simulation (LES) for the simulation 
of high-speed reacting flows [16]. LES methods calculate the large scale features of the flowfield 
and require modeling only to describe the small scale features of the flow utilizing a subgrid scale 
model. Subgrid scale models using filtered density functions (FDF) have shown great promise in 
improving the state-of-the-art in combustor modeling [17]. The FDF method is the counterpart of 
the probability density function (PDF) methodology, which has been long-established in RANS 
calculations of turbulent combustion [18]. The axisymmetric version of the SPARK RANS code 
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Figure 10: Comparison of RELIEF velocity data with simulation results 



was used to simulate the flowfield in the helium/oxygen center nozzle and the outer air nozzle of 
the coaxial jet mixing experiment described in the previous experiment. Details of the code are 
given in reference [19]. The analysis of the experiment was begun by first solving for the flowfield 
in the center and outer nozzles. Results obtained at the end of each nozzle were then used to spec- 
ify the supersonic inflow conditions for the downstream domain beyond the nozzles where mixing 
of the jets occurred. Details regarding the nozzle calculations and the mixing region downstream 
of the nozzles are given in reference [20] . 

A comparison of the measured helium mass fraction data with the simulation results at several 
stations downstream of the nozzles is given in Figure 8. Agreement between the simulation and 
the data is very good at each station. The code somewhat overpredicts the mixing near the center- 
line at the x = 0.12 m station, although the prediction improves with increasing radial distance. A 
comparison of measured pitot pressures with the simulation is shown in Figure 9. Agreement is 
good in the region of the air coflowing jet, but the simulation somewhat overpredicts the pitot pres- 
sure in the helium-air mixing region. The comparison with the experimental data differs at large 
radial distances greater than 0.025 m as the code does not consider the effects of the laboratory air 
entrained by the coaxial air jet. The RELIEF streamwise velocity data is compared with the simu- 
lation in Figure 10. The prediction agrees well with the data at the first three stations and slightly 
overpredicts the data at the remaining stations near the centerline. The simulation somewhat un- 
derpredicts the the velocity at the final three stations in the mixing region between the helium and 
air coflowing jets in agreement with the pitot pressure results. 


5 Concluding Remarks 

This paper has described work to develop phenomenological models for Reynolds averaged Navier- 
Stokes codes, subgrid scale models for use in large-eddy simulation and reduced-kinetics mecha- 
nisms to model hydrogen-air and hydrocarbon-air reaction in propulsion applications. Fundamen- 
tal experiments are being performed to provide data that will be used in the development and re- 
finement of these models. Experimental data is extracted from the experiments using nonintrusive 
diagnostics that allow accurate simultaneous measurement of temperature, species, and up to three 
components of velocity in supersonic flow without changing the character of the flowfield. Once 
the databases are in hand, the data will be analyzed using a response surface methodology [21-23] 
that provides an efficient means of determining critical parameters in a chosen model. The models 
will then be incorporated into combustion codes used in engine flowpath analysis and design. 
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